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We study the real-space entanglement spectrum for fractional quantum Hall systems, which main- 
tains locality along the spatial cut, and provide evidence that it possesses a scaling property. We 
also consider the closely-related particle entanglement spectrum, and carry out the Schmidt decom- 
position of the Laughlin state analytically at large size. 

PACS numbers: 



It is by now widely accepted that topological phases 
of matter Q defy the paradigmatic classification based 
on local order parameters and broken symmetries. In- 
stead, quantum information concepts, such as quantum 
entanglement, have yielded valuable insight on this topic 
in recent years, and are believed to be good probes of 
topological order. For instance, the topological entangle- 
ment entropy @, 0] obtained from a spatial bipartition 
of the ground state of a fully-gapped Hamiltonian in two 
dimensions is one characteristic of a topological phase. 
At the end of Ref. 0, Kitaev and Preskill (KP) showed 
that the same result can be obtained by assuming that 
the reduced density matrix has the form of the thermal 
density matrix of a gapless chiral energy spectrum, as 
would occur at an edge of a quantum Hall (QH) system. 

To be precise, a bipartition of a quantum system is 
defined when the Hilbert space factors into two parts, 
H = Ha ®Hb- Then the (normalized) ground state has 
a Schmidt decomposition [4|: 

\^)=Y j e-^ 2 \i) A , i )®\i> Bii ), (1) 
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where the e~^/ 2 's (e - ^/ 2 > 0) are the Schmidt singular 
values, and \ipA,i) (iV'-B.i)) are an orthonormal set in Ha 
(Hb)- Equivalently, the reduced density matrix pA = 
tr-H B IV 7 ) (V'l nas eigenvalues e~^ i , and J2i = (V'lV') 
1 . The set of & is called the entanglement spectrum (ES). 

For extended homogeneous systems defined as many- 
particle systems, lattice models, or field theories, with 
fairly short-range interactions, a natural way to partition 
is with a cut in position space (in continuous systems), 
or between a row or plane of sites (in a lattice model). 
It is for such types of partition that the results of Refs. 
0, 13] were derived. But studies of fractional QH states, 
in which the particles are confined to the lowest Landau 
level (LLL), have usually used other bipartitions. One of 
these is the orbital partition (OP) [f|, in which a basis 
set for the LLL single-particle states is divided into two 
subsets, and this leads to a bipartition of the iV-particle 
Hilbert space also. If the basis set has the form that is 
natural in the Landau gauge (we will describe it for the 
plane here for simplicity, though later we use a compact 
geometry), then basis states are eigenstates of transla- 



tions in the y direction (with momentum eigenvalue k y ), 
and a Gaussian form in the x coordinate (with center at 
k y ), so that partitioning into, say, the subsets k y < and 
k y > corresponds roughly to a spatial cut along the y 
axis. For the LLL on a sphere, the corresponding scheme 
uses eigenstates of the 3-component of angular momen- 
tum, and corresponds roughly to a spatial cut along a 
line of latitude. 

Li and Haldane (LH) @ studied the ES of QH states 
using OP on the sphere. They argued that the ES con- 
tains a low-lying part in which the multiplicities are re- 
lated in a universal way to those in the chiral confor- 
mal field theory (CFT) which describes the edge physics 
of the topological phase. This low-lying part is usually 
well-separated from the rest of the ES, with a gap going 
to infinity for certain model wave functions such as the 
Laughlin state 0] ■ LH proposed to use the low- lying part 
as a diagnostic tool for numerical calculations. 

In this paper, we examine the direct real-space par- 
tition (RSP) in QH systems, defined in coordinate (not 
momentum) space, just as it is for non-QH systems. This 
bipartition by definition maintains locality along the cut 
in, say, the y direction — that is, correlations (described 
by p A ) along the cut are short range. Such locality, which 
(as we will show) does not hold for the OP that involves 
instead a bound on k y , is extremely helpful in theoretical 
analysis. First, for RSP of Laughlin states, we provide 
numerical evidence that a "scaling form" motivated by 
locality holds for the ES. Stated loosely, this says that 
as the thermodynamic limit is approached, not only the 
multiplicities but also the eigenvalues & in the low-lying 
ES approach the energy spectrum of a local field theory 
on the cut, which is similar to the edge theory (for the 
trial states we study, this is a chiral CFT) . This stronger 
form of the LH proposal resembles the idea of KP. It is 
not clear if the topological entanglement entropy 0, Hf is 
obtained if this scaling form does not hold. Second, we 
also consider another bipartition for QH systems, called 
particle partition (PP) [5j, for which we extend the orig- 
inal definition, and show its close relation with RSP. Fi- 
nally, we show for RSP and PP that the multiplicities for 
the Laughlin states obey the scaling form, and find the 
Schmidt decomposition explicitly for PP. 
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Long-range correlations along the cut in OP — We 
use the example of the equal-time Green's function (off- 
diagonal one-particle density matrix). In a large system 
of particles in the LLL in the x-y plane, this takes a form 
independent of the details of the ground state other than 
the filling factor v, provided the state is translationally 
and rotationally invariant: 



(iP\ft(x',y')i>(x,yM = _ e (^') 2 /4-, 2 /2-, 
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(2) 



in the Landau gauge, where we use z — x + iy, z' = 
x' + iy' , and ip(x,y) is the destruction field operator at 
z (in this paragraph only, we set the magnetic length to 
1). This function falls off as a Gaussian in \z — z'\ in all 
directions. It can be derived by expanding ip(x, y) in the 
(continuous) basis of LLL momentum eigenstates in the 
y direction, which are oc e~^ x ~ ky ^ / 2 + lk vV : 

(^(x',y')i,(x,yM 

dk y e k ^ z+ ~">- k2 y~ x2 l 2 - x ' 2 > 2 . (3) 
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For OP, we must restrict the integral to k y in part A, say 
k y < 0. Due to the discontinuity in the integrand of the 
Fourier transform, we find 
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k v {z+z')-k 2 -x 2 /2-x' 2 /2 



2 /2-x' 2 /2 



2w 3 / 2 i(y - y') 
(4) 

as \y — y'\ — > 00; the long-range part has largest ampli- 
tude at x = x 1 — 0, the "cut" . We cannot evaluate other 
correlations as easily, but they will have related asymp- 
totic forms. Thus in OP, correlations evaluated within 
part A are long-ranged in coordinate space. 

Quantum Hall states on the sphere and RSP — The 
single-particle states in the LLL on the sphere S 2 pierced 
by a rotationally-invariant magnetic field with a total of 
N$ flux quanta (N^ > an integer) [8[ can be described 
using stereographic projection to the plane. An orthonor- 



mal basis set of LLL states ip r , 



(m 



0, 1, 



N^) with 



3-component of angular momentum L z = /2 — m on 



the sphere is thus given by 



r m V(z)/2 



up to a 



normalization factor (defined using the L norm on the 



x + iy plane), where the function V(z) 
<y(z) 1 



(1 + M 2 ) 



2+Na 



(5) 



arises from the coordinate transformation from the 
sphere. The northern and southern hemispheres corre- 
spond to \z\ < 1 and \z\ > 1 respectively. 

For a system composed of many identical particles, a 
bipartition can be defined by first dividing the single par- 
ticle Hilbert space Hi into two subspaces as a direct sum: 
"Hi = %xa®'Hib- This then induces a corresponding bi- 
partition of the TV-particle space, in the form 

JV 

®T-Ln b m, (6) 
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FIG. 1: (color online) Low-lying RSP ES of the v = 1 filled 
LLL, for N — 30 particles, and Na = 15. Lower inset: lowest 
levels. Upper inset: low-lying ES versus Na for N = 16. 



where Na (Nb) is the number of particles in part A (B). 
For RSP, "Hi is an L 2 space of all normalizable func- 
tions over a manifold (the sphere in our case), and the 
two parts are simply the subspaces of normalizable func- 
tions supported on a subregion A or on its complement 
B. Then any single-particle wavefunction can be decom- 
posed into two orthogonal parts, each part lying in one 
subspace: ip(r) = tp A (r) + ip B (r). The corresponding 
bipartition of the iV-body Hilbert space is what we call 
RSP. In this paper, we usually take region A to be the 
northern hemisphere, and B the southern. We empha- 
size that for the QH problem, the L 2 space in which we 
decompose contains all the Landau levels and is infinite 
dimensional. Even though we consider ground states in 
the LLL, Ha (and Hb) in the RSP decomposition include 
states that are not LLL states. The Schmidt decomposi- 
tion is nonetheless of finite rank for a finite-size system. 

RSP for the v = 1 filled LLL — This case can be solved 
analytically as follows Q. Let c m be fermion operators 
which destroy particles in the LLL basis states functions 

ip m . Then the filled LLL is IImio4J )> where 1°) is 
the vacuum with no particles. For the Schmidt decom- 
position, we define operators and which destroy 
particles in the normalized states obtained from ip m by 
multiplying it by 0(1 — |z|) or @{\z\ — 1) respectively, and 
normalizing. Then we have 



Cm — ^m^jn ~\~ /^mC m , (T) 

where a m > 0, (3 m > are the restricted norms: 



= J d 2 z|^ m | 2 0(l 



1*1), 



(8) 



and similarly for j3 m using region B. Substituting this 
into the ground state yields the Schmidt decomposition 
immediately. For each occupied basis state, the particle 
can be in either part A or part B, and these occur with 
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FIG. 2: (color online) Similar to Fig. [TJ but for the v = 1/3 
Laughlin state with N = 12, Na = 6, and upper inset shows 
minimum £ only. 

weights a m or (3 m . Hence each singular value e - ^/ 2 is 
given by a product of a m 's or /3 m 's, one for each m. It 
follows that the terms i in the Schmidt decomposition 
can be labeled by a set of "occupation numbers" n m = 
or 1 for each m, and the many-particle pseudoenergies & 
are given by a constant plus a sum of single-particle pseu- 
doenergies for the occupied modes: J2m=o n m£m where 

e m = -ln(a 2 m /pl). (9) 

e m is monotonically increasing with m, so the lowest £ is 
found by occupying all m with m < and leaving 

the rest unoccupied (there is a zero mode Sns/2 = for 
iV^ even). The ES thus has exactly the form of the edge 
of the v = 1 state. Here JVa = £ m n ™> -/Vb = N — N A . 
A sample spectrum is plotted in Fig. [TJ for N A = N/2 
as a function of Lf, the angular momentum of part A. 
The linear dependence containing a velocity, which comes 
from the linearity of e m near m = 0, is evident. The 
full ES versus Na is also shown. The lower edge has a 
parabolic form, which we discuss below. 

In performing the decomposition, we must reorder 
fermion operators, so that those referring to part B are 
to the right of those for part A, because the tensor prod- 
uct in eq. (J6j> is Z2 graded. This introduces some minus 
signs, which become important for correlated states. 

For the OP, the ES of the filled LLL contains only a 
single state, with Na and determined by filling all 
orbitals in part A. Thus the LH conjecture about the 
multiplicities does not hold in this case. 

Scaling conjecture — Turning now to fractional QH 
states, we can state what we expect for the behavior 
of the ES in RSP. In the low-lying branch for generic 
ground states, or in the full ES for trial ground states, 
there is usually a single non-degenerate lowest pseudoen- 
ergy level £0 at some values Nao and Lf of the good 
quantum numbers Na and Lf; these values depend on 



the system size N. Let us define quantum numbers and 
pseudoenergies AN a, ALf and A£j by subtracting off 
these values. We will say the ES has the asymptotic scal- 
ing property if for all AN a and ALf, as N — > 00, the 
set of A£i approach the energy levels (minus the lowest 
energy) of a Hamiltonian that is the integral of a sum 
(with size-independent coefficients) of local operators in 
a 1 + 1-dimensional field theory on a circle. In particular, 
for certain trial states, such as the Laughlin and Moore- 
Read states, the field theory will be a chiral CFT that is 
the same one as for the edge theory. In the simplest case, 
the leading local operator allowed by symmetry will be 
the stress tensor T(z), and then the low-lying spectrum 
will collapse onto a single straight line versus ALf for 
each AN a- The coefficient is a "velocity" v that actually 
has dimensions of length; if it is non-zero we have 

A& - vLo/R (10) 

where Lq = — AL^ + 0{AN a) is the ith eigenvalue of the 
zcroth Virasoro generator, and R is the radius of part A; 
in our scaling v/R = 0(N A 1 ). We emphasize that this 
scaling property encompasses both (i) the multiplicities 
and (ii) the pseudoenergies at each AN a and ALf. 
LH's conjecture was that scaling holds for the multiplic- 
ities only. We conjecture that the scaling property holds 
in full for RSP ES of a gapped topological phase. 

Numerical results — General trial states for a partially- 
filled LLL are linear combinations of Slater determinants, 
and the calculations resemble the steps above for the 
filled LLL, but the final Schmidt decomposition was per- 
formed numerically. Computations for RSP are much 
more costly than for OP of the same state. We present 
results for the Laughlin v = 1/3 state in Fig. [21 The 
ES is strikingly linear at small AL^, and multiplicities 
there are the same as in the edge of an infinite system for 
the lowest values, and the same as for OP over a larger 
range. We also show the lowest pseudoenergies as a func- 
tion of AiV^. The parabolic form should be compared 
with what is expected for the CFT of the edge theory (a 
chiral compactified scalar field [l(| ) , which by the scaling 
conjecture should be 

At (AM , v Q(ANa) 2 

for filling factor v = 1/Q (including v = 1). For the data 
in Fig. [51 this works quite well, by taking v /R from the 
spectrum at fixed AN a- It gives a numerical estimate for 
the "compactification radius" for the scalar field, which 
should equal Q. Moreover, the general form of the spec- 
trum is similar to that for v — 1, for which we know ana- 
lytically that the levels eventually collapse onto a straight 
line. While finite-size limitations have prevented us from 
making scaling plots to study the approach to scaling, it 
is at least plausible that it will hold based on the trends 
visible at these and smaller sizes. 
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FIG. 3: (color online) Same as Fig. [2] but for the Coulomb 
interaction. 



FIG. 4: (color online) Low- lying ES for the PP in the v = 1/3 
Laughlin state, for N — 12 particles, including all Na values, 
versus L A , not Lf . 



For the LLL Coulomb interaction, the low-lying ES is 
similar to that of the Laughlin state — see Fig. [3] We 
have also examined the ES for RSPs of the Laughlin 
states with cuts along different lines of latitude (not 
shown); these are always similar, but the minimum of 
the parabola of A£ m ; n versus ANa is shifted. 

Particle Partition — Another bipartition is the particle 
partition (PP) 5] in which one divides the TV particles 
into two sets of fixed sizes Na and Nb, and calculates 
the entanglement. The ES in PP has been considered 
recently 11] . We will first extend the definition so that 
it yields all values of Na, as for OP and RSP. We intro- 
duce fictitious "pseudospin" variables for each particle; 
the two orthogonal pseudospin states will be called \A) 
and \B). To map a spinless ground state into the larger 
space that includes pseudospin, we assign to each particle 
j the pseudospin state (\A)j + \B)j)/y/2 (thus the state 
is still totally antisymmetric). We define a bipartition 
of this state in which the particles with pseudospin \A) 
constitute part A, and those with \B), part B. If oper- 
ators c" n destroy particles in orbital m and pseudospin 
a = A, B, then the decomposition ([7} again applies (c m 
being for the superposition pseudospin state), but now 
with a m = /3 m = 1/V2 for all m. This is our extended 
PP; the part with any fixed Na is essentially the old PP. 
This scheme can be generalized by utilizing arbitrary val- 
ues of a m , (3 m (Q!to + Pm = 1); OP also fits into this by 
using values a m = or 1. Now changing the set of a m to 
another set of values (say, 1 / \/2 for all m) defines a linear 
map Hn a ,a ~ > Wn a ,a (more precisely, of the subspaces 
spanned by the basis states), which is invertible if all a m 
are non-zero (and similarly for part B); it follows that 
the Schmidt rank in RSP is always the same as in PP, 
for each Na and L A (but generally larger than in OP). 

For PP, rotational invariance is preserved, so we may 
use the total angular momentum L A of part A to label 
multiplets, each of which has degeneracy 2L A + 1 — see 



Fig. 2J Unlike RSP, if £ m i n is plotted versus Na, it yields 
an inverted parabola. For PP of the filled LLL (or for 
any single Slater determinant), — Nln2 for all i. 

Explicit Schmidt decomposition of Laughlin states — 
We perform the PP of the Laughlin state ^>q [with 
Ncf, — Q(N — 1)] by letting the coordinates for part A 
be zt, ... , zn a , and for part B, w\, . . . , Wn b . Then the 
wavefunction is 



i,l i<j k<l 



x e z 



oc exp 



] 3— m3 — m J 

m>0 / 



*q({z})*q(M),(12) 



where we have used the following notations: 

l<i<j<N A 



k=l l<i<j<N B 
N A 




and note (z — w)® — {— w)® exp(— Q J2m=i z ™ w m /m) 
for \z\ < \w\. Note that ^q({z}) is a droplet centered 
on the north pole, and ^q({w}) is centered on the south 
pole. Each factor j A m or j^ m 's creates a density excita- 
tion on the corresponding edge [hJ. The terms that are 
obtained by expanding the exponential in (1121) are of the 
form 



II 



(13) 



m>0 
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times similar for B. A linearly independent set for part A 
is those with n m = for m > Na and J2 m mn m < QNb- 
Hence, this proves easily that in the limit Na, Nb —> oo, 
the count of levels for each ALf and ANa (defined using 
Nao = N/2, Lf = maxLf at Nao) matches that of the 
edge theory. 

To complete the Schmidt decomposition, we need to 
know if these basis sets for parts A and B are orthonor- 
mal in the limit. This follows from existing results, in 
Ref. to leading order, and Ref. [IH to higher order. 
This implies that for PP, the in the scaling region 
obey the scaling property, with vanishing coefficient in 

— 1/2 

the 0(N A ) term, that is vanishing velocity, for all Q 
(this can also be shown for AN a ^ 0). For Q = 1, we 
know this holds to all orders. For all Q, if the scaling 
property holds, then the vanishing of the velocity follows 
from rotational invariance of PP. For RSP, we may imag- 
ine that a non-zero velocity appears by deforming the PP 
results due to the restriction on particle coordinates when 
calculatin g th e overlaps. We present a detailed analysis 
elsewhere 13 1. 

Conclusion — The RSP leads to very attractive results 
for the ES, and appears to obey the scaling property 
that results from locality. The PP is closely related, and 
allows an analytical treatment for the Laughlin state, in 
which the limiting ES in the scaling region has vanishing 
velocity. 
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